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Optically mosaic crystals of zinc have been grown by the Czochralski-Gomperz method. 
Three types have been distinguished. During growth mosaic crystals frequently change from 
one type to another. The conditions of growth under which they were obtained are described 
briefly. As yet their genesis and their evolutions during growth cannot be controlled or pre- 


determined. 


URING experiments on the growth of single 

crystals of zinc by the Czochralski- 
Gomperz method the writer has obtained several 
hundred specimens with “optically mosaic’”’ 
structure.! The term ‘‘mosaic’’ may properly be 
applied, as it is in the literature of nonmetallic 
crystal physics,” to a crystal having a ‘“‘broken”’ 
cleavage surface which consists of plane areas 
each slightly differently oriented with respect to 
its neighbors. When these elementary areas are 
large enough to be observed with a low power 
microscope or the unaided eye, the crystal may 
be called optically mosaic. 

The three distinct types of. such mosaics en- 
countered are illustrated by Figs. 1 to 3, which 
are photomicrographs of unetched, unpolished, 
natural cleavage surfaces.* In type A, Fig. 1, the 
components of the mosaic pattern have irregular 
boundaries and are small compared to the whole 
surface. An entire cleavage surface may consist 


*The writer has just learned that similar mosaic zinc 
crystals were grown by Dr. James E. Brock about three 
years ago in the Physical Laboratory at Purdue University. 
The method of growth, however, was not the Czochralski- 
Gomperz method but the one devised by Bridgman. 

* Wien-Harms, Ott, Handb. d. Exp. Physik 7, 2nd part, 
101, 102. 

* Basal planes of the zinc crystal. 


of mosaic elements, as in the figure, or the mosaic 
pattern may be confined to only a part of it. 
The differences in orientation, shown in the 
photograph by the differences in brightness of 
neighboring areas, are usually less than a degree. 
That each area is itself plane is shown by the 
fact that any given one always appears equally 
bright all over. In type B, Fig. 2, there are large 
areas bounded usually by straight lines and occa- 
sionally by smooth curves. The angle of tilt be- 
tween two neighboring areas may amount to 1 or 
2 degrees. In type C, Fig. 3, the elementary 
mosaic areas are bounded by nearly parallel lines 
which to the unaided eye appear straight, but 
which under greater magnification are seen to be 
jagged. For crystals having orientations of 35° or 
more these bounding lines are always very nearly 
parallel to the major axis of the elliptical cleavage 
surfaces of the writer’s circular cross section 
crystals. 

A single mosaic type does not necessarily per- 
sist throughout a crystal of 8-10 cm length. Fre- 
quently the ends of a specimen may be of differ- 
ent types with a region of transition between, or 
the same type, neither mosaic or single, may oc- 
cur at both ends, the middle being quite differ- 
ent. Changes from one type to another occur 
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Fic. 1. Type A mosaic crystal. 
Diameter of area photographed = 
D=3.3 mm. 


apparently at any speed, sometimes taking place 
during only a few mm of growth, and sometimes 
requiring 10 or 20 cm of growth for their com- 
pletion. In the majority of cases the course of 
such changes is quite unpredictable. However, 
one typical evolution, a case of change from single 
to type C structure and back to single, has been 
observed. Its various stages, which are quite 
invariable in their sequence, are sketched in 
Figs. 4 and 5. The positions along the crystal of 
successive cleavages (for a crystal having an 
orientation of approximately 45°) and the ap- 
pearances of corresponding surfaces are repre- 
sented by a, b, c, etc., of the two figures. Section 
a (near the top or nucleus end of the crystal) is 
optically single. The next section, b, shows a 
small area of slightly different orientation mak- 
ing its appearance at the upper end of the major 
axis b’. In the following sections this area in- 
creases in size, growing first mainly along the 
major and axis, then becoming wider until there 
are three distinct parallel-edged areas. Two of 
these are then crowded out by degrees until 
_ finally we again have an apparently good surface 
(h) which may have the same orientation as (a) 
or one slightly different. While scores of such evo- 
lutions have been observed, not a single one has 
been found in which the invading area made its 
appearance at the lower end of the major axis 
(at b in Fig. 5). Nor does it ever make its ap- 
pearance on the sides of the elliptical periphery, 
though it may grow down one side of the major 
axis turning the crystal into a double of type C, 
as illustrated by Fig. 6, a to e, instead of the 
triple as in the previous one. 
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Fic. 2. Type B. D=6.4 mm. 


Fic. 3. Type C. D=4.4 mm. 


The apparatus was essentially that of Hoyem‘ 
except that the crystals were grown in air. The © 
growing procedure was at first as near as possible _ 
that of Hoyem and Tyndall,’ though later modi- 
fications were made which need not be described 
here. Various speeds between 0.1 cm/min. and 
8 cm/min. were used. Great care was taken to 
ensure uniform motion of the drawing apparatus 
and to prevent vibration. The specimens were _ 
from 8 to 30 cm in length and from 0.2 cm to 1.1 
cm in diameter. The diameter could with care be 
kept uniform to within 0.01 cm. Mosaics have 
been grown over a temperature range extending 
50° above the melting point of zinc and over the 
whole orientation range. During growth the 
temperature of the molten zinc in the crucible 
immediately below the growing crystal was main- © 
tained constant to within 0.5°C. The rate of dis- 
charge of the air blast directed from a ring nozzle 
to the crystal at a point 1.5 cm above the melt 
was constant. 

Nine different brands or lots of zinc were used. 
With one exception they were all of reagent grade 
and some certainly contained not more than 0.001 
percent of impurity. Mosaics were obtained from 
all these zincs although there were some slight — 
differences in the type produced, etc., in the — 
various lots. 

Mosaics could appear when either a single or a 
mosaic nucleus was used to inoculate the growing . 
crystal. In spite of the maintenance of very rigid © 
control of all the factors ordinarily considered im- © 
portant in crystal growth, it has thus far been 


4 Hoyem, Phys. Rev. 38, 1358 (1931). 
5 Hoyem and Tyndall, Phys. Rev. 33, 81 (1929). 
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Fig. 4 


OF ZINC 3 


Fig.6 


Fics. 4-6. Stages of typical evolutions from single to type C structure and back to single. 


impossible to control the appearance of mosaics 
or prevent the occurrence of evolutions during 
growth. Both mosaic and single crystals have been 
obtained for most lots of zinc under apparently 
identical conditions of growth. The investigation is 
being continued in the hope of discovering why 
the writer’s experience is so different from that 


of previous workers, in this laboratory and else- 
where, who grew single crystals with no difficulty 
and encountered no mosaics. 

This work is being done under the direction of 
Professor E. P. T. Tyndall, to whom the writer 
is greatly indebted for many suggestions and for 
encouragement. 
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Diffusion Through a Membrane 


CoLin BARNES, McLennan Laboratory, University of Toronto 
(Received July 21, 1933) 


An exact solution is obtained for the problem of diffusion 
through a membrane separating two stirred solutions of 
different concentrations. For a diffusion cell in which the 
two liquid containers have the same volume and for which 
the ratio \ (“‘effective’’ volume of membrane : volume of 
liquid) is small compared with unity, the solution gives 
explicit expressions for the concentrations at any time in 
terms of the initial distribution of concentration in the 
membrane. Two cases are treated: (1) zero initial concen- 
tration in the membrane and in one of the liquids; (2) an 
initial linear concentration gradient connecting the concen- 


trations in the two containers. The usual assumption of a 
linear gradient across the membrane during the whole 
course of diffusion will not lead to an appreciable error in 
a determination of a diffusion constant unless \ exceeds 
0.1 or the time of diffusion is much shorter than in the 
usual practice. Since \ is about 0.02 in the usual type of 
diffusion cell the use of the simple expression based on a 
linear gradient cannot be in error. If, however, the pre- 
liminary treatment of the membrane fails to yield an 
initial linear gradient in it, an appreciable error will be 
made even with X less than 0.02. 


1. INTRODUCTION . 


SIMPLE and precise method for the 
determination of diffusion constants in- 
volves the use of a diffusion cell in which a 
porous membrane separates two homogeneous 
liquids.! Thus a solution of concentration 7 is 
placed on one side A of the membrane and pure 
solvent on the other side B. The membrane is 
given a preliminary treatment to justify some 
statement about the initial distribution of the 
solute in it. Diffusion is allowed to proceed for 
a definite time, the liquids being kept homo- 
geneous by (say) stirring. The concentrations in 
A and B change from their original values 7%, 
0 to different values v, w in the time ¢, and from 
a knowledge of these final values it is desired to 
find the diffusion constant of the solute. 
It is usual to assume that the preliminary 
treatment of the membrane can be made to give 
an initial linear distribution across it and that 


linearity is maintained during diffusion.” If the 
two containers A, B have the same volume D 
and the membrane has an effective volume a 
(determined by the distribution and structure 
of the pores), the assumptions give at once the 
relations 


where is the diffusion constant and 
the amount of substance in the membrane 
remains constant. The cell-constant d/a’® is 
usually determined by calibrating the cell with 
a solute of known diffusion constant. 

Whatever may be the opinion regarding the 
possibility of producing an initial linear gradient 
across the membrane, the assumption of a linear 
gradient during diffusion contradicts the exact 
equations for the problem. An exact solution 
from given initial conditions may, however, be 
found.* 


2. GENERAL SOLUTION OF THE PROBLEM 


Let the membrane be bounded by the planes x =0, x =a. Well-stirred solutions of initial concen- 
trations v, wo are in contact with the two sides of the membrane; the first solution extends from 
x=0 to x= —b and the second from x =a to x=a-+tc. The initial distribution of concentration across 


1 J. H. Northrop and M. L. Anson, J. Gen. Physiol. 12, 
543 (1929). 

2J. W. McBain and T. H. Liu, J. Am. Chem. Soc. 53, 
59 (1931). 


3H. W. March and W. Weaver, Phys. Rev. 31, 1072 
(1928), have solved a related problem of the diffusion from 
a gel into a liquid. 
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the membrane is mo(x). After a time ¢ the concentrations in the two solutions will be denoted by 
v(t), w(t) and the distribution in the membrane by u(x, ¢). 
The equation 
dt=kd?u/ dx’, (2) 


where & is the diffusion constant of the solute, is to be solved subject to the following initial and 
boundary conditions: 


u(x, t)—>uo(x), t0, (3) 
u(0, t) =v(t), u(a, t)=wi(t), t>0; 7(0) =v, w(0) (4) 
dv/dt = dw /dt= t>0, (5) 


the Eqs. (5) expressing the proportionality of rate of increase of concentration of solution to rate of 
supply of solute from the membrane. 


The solution of (2) which vanishes at x =0 and at x=a and which satisfies (3) is 


= sin (n2x/a), (6) 
5 
where 
Cu=(2/a) sin (ney /a)dy (7) 
0 
and 
an =kn?x?/a?. (8) 


The solution of (2) which corresponds to zero initial concentration in the membrane and the faces 
maintained, respectively, in contact with solutions of concentrations v(t), w(t) is known to be‘ 


(u(r) (0/ dt) Fy(x, t— 7) + w(7)(0/dt) Fo(x, t—7)} dr, (9) 
where 
t) =1—x/a—(2/7) sin (u7x/a) (10) 
and 
F3(x, t) =x/a+(2/7) 1)"/n)e—*! sin (n2x/a). (11) 


We may write 
t 
-{ Ar) Fi(x, t— 7) +w(7)(0/d7) Fo(x, t— 7) 
0 
t 
=voF\(x, t)+woFo(x, t) +f (7) Fi(x, t— 7) (7) Fe(x, t— 17) }dr, (12) 
0 


and the solution of (2) which satisfies (3) and (4) is simply u=u’+w’’. The value of u so found from 
(6) and (12) gives 


—a(du/dx) =v—w+2> cos (nrx/a)+2>¥ cos a) enone —(—1)"w’)dr, 


where 
(13) 


Eqs. (5) now give, respectively, 


* See, e.g., Carslaw, Mathematical Theory of Heat Conduction in Solids, 1921, pp. 18, 68. 
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and 


These equations will be satisfied identically by putting® 


dv/dt= > ***, dw/dt=> (14) 
if the A;, B;, &; satisfy the equations 
A nl /(an—&)) +220 n(—1)"/(an—&i)), (15) 
Dn= (17) 


Eqs. (14) will provide the solution of the problem if (15), (16), (17) lead to unique values for the 
numbers A;, B;, With 


(18) 
the last three equations become 
A,(1/i—cot i) = —B,(1/2; sin Zi), (15’) 
B1/dA2—cot 2;/2;) = —A,(1/2; sin Zi), (16’) 
D,! =kD,/a@ => (Ai—(—1)"B,)/(n? —32?), (17’) 
with 
A, =a/b, ho=a/c. (19) 


Eliminating A;, B; from (15’), (16’) we get the equation 
— (Ai +A2)z cot Z—A,A2=0, (20) 


the real positive roots of which give the numbers §;. 

From (15’) or (16’) B; is given in terms of A;, &;, and the infinite set of linear equations (17’) 
provide the numbers A;. Regarding the A;, B;, &; as known, the concentrations of the two solutions 
are 


(21) 
where @, the final uniform concentration, is given by 
Uo(x)dx+ woe. (22) 
0 
From (21), 
>A i/ > =u —wo; (23) 


one of these is needed for the numerical evaluation of the A;, Bi. 
By means of (21) for v, w, the concentration distribution in the membrane is eiiail from (6), 


(12) to be i 
u(x, sin (Ai—(— (24) 


It may be verified that this expression satisfies (2)—(5) when due account is taken of (21) and the 
conditions involved in (15’)—(17’). The arguments necessary to make the above proof rigorous are 
elementary and readily supplied; they are omitted as providing nothing of interest. 


5 March and Weaver, reference 3; E. T. Whittaker, Proc. Roy. Soc. (London) A94, 367 (1917). 


7 
; 
: | 
|| 
| 
| 


t) 


DIFFUSION THROUGH A MEMBRANE 7 


3. APPLICATION TO DIFFUSION EXPERIMENTS 


In a usual type of apparatus constructed for measuring k the quantities \;=a b and \»=a_¢ are 
equal and small compared with unity (of the order of 1/50). The general solution may be adapted 
at once to this case and all the accuracy desirable in practical cases is obtained by disregarding powers 
of \(=A:=As) greater than the second, as negligible compared with unity. An equality sign will 
frequently be used where such an approximation has been made. 

Let the positive roots of (20) with 4; =.= be 2;, 7=0, 1, ---, arranged in order of magnitude; 
standard processes® for finding the roots of such an equation give 


So = (25) 
+---, i=1,2,---. (26) 


When A, =e, (15’), (16’) show that 


and the linear equations for the A; become 
(nx? —3?). (27) 


The values of the A; depend on the initial concentration in the membrane, i.e., on u(x) of (3). 


Case 1 
A case of practical interest is that in which u(x) =0, wy) =0; D,,’ = kvo/a®. Here the Eqs. (27) become 
= aniA i 
1 


with 
= kvo/a?—[(1+(—1)") / (n?2? — 20?) JA, / —2?). 
The solution is’ 
2?) Ao a?) 27), for 7 even, 
A;= a?) (1—2d/7? 2’), for 7 odd. 


With the roots (26) these values of the A; give 
A, 52 =(4/360)A 492/360), 
1 


and from the first of Eqs. (23), (25) and %—v = —(v/2)(1+A/2—?/4), we get 
(Ay 2d)(1+A/6+ 22/90) + (ke /2a2)(1 —¥/6+ 22/45) =0. 


Hence 
Ao/ = — 
A & = — (vo, 2)(40/P 2”) (1 2”) i even or odd. 
Thus 
Vo A o 4A OX 
2 2 4 6 60 1 
Vo A » OX 
(1 —-+— — (1 --+— / 6+? / 45) —— 
2 2 4 6 60 1 Hy? 1? 3 


° See, e.g., Geiger-Scheel, Handb. d. Physik 3, 600 (1928). 
7 Only two terms in the solution-series are needed: a@;;A;=E;—L,a;,C;/a:,. For existence theorems see H. v. Koch, 
Jahresb. d. d. Math.-Verein. 22, 289 (1913). oo 
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Case 2 


Another case of interest is given by uo(x) =1(1—x/a), =0; C,(Eq. (6)) = D,’ =0. Here 
the concentration gradient in the membrane is initially linear. In the notation of Case 1, 


E,=((1+(—1)")/ x? — 20?) 


and the a,; are as before. The solution for the A; is® 


A;=(4X for 7 even, 
A;=0, for 7 odd. 
With these values the first of (23) gives, with 7 =7»/2, 


ki 2a*+ (Ay 6+. 90) =0. 
Hence 


Ay b= —(vo/2)(1— 180); Ay &=---=0; As; £o;= — (vo (1/(27)4). 
The concentrations v, w are thus 


2 18( ) 


1 


2 180 2rii 


13 


These equations indicate the error involved in the use of (1). The change from the exponent 
a’ in (1) to the exponent A(1—A_ 64+ 45), a? is of no experimental significance; either quantity 
represents a cell-constant determined by calibration. The error in (1) enters through the neglect 
of \*/180 in comparison with unity and the neglect of the terms following the first exponential. 
In the usual diffusion cell \ is of the order of 0.02 and 2dkt, a? is of the order 1 (diffusion proceeds 
until v is about 37). The terms following the first exponential are entirely negligible and the error 
introduced by the term \?/180 corresponds to an error of less than 1 in 10° in the concentration 
change. The deviations from (1) do not become appreciable until \ exceeds 0.1 and the time of 
diffusion is taken so short that the remaining terms of the series must be taken into account. Hence 
with the diffusion cell as ordinarily used no appreciable error will be introduced by the use of the 
simple expressions (1). A much greater error is to be expected if the initial gradient is not linear 
(Case 1). 

I wish to thank Professor A. R. Gordon of the Department of Chemistry, University of Toronto, 
for suggesting the utility of an investigation of this problem. 


5 See footnote 7; only one term in the solution-series is now needed. 


‘ 
. 
$4 
te 
= 
re 
| 


JANUARY, 


1934 


PHYSICS 


VOLUME 5 


A Thermoelectric Metal-Crystal Analyzer 


M. F. HASLER, California Institute of Technology, Pasadena, California 


A device is described which allows the measurement of the thermoelectric ¢.m.f. at room 
temperature for each increment of length of a metal test bar. A photographic arrangement 
permits a complete thermoelectric record of such a bar to be produced. 


INTRODUCTION 


T is a well-known fact that the thermoelectric 

e.m.f. is a very sensitive function of the 
physical state of the sample upon which it is 
measured. Mechanical work within the elastic 
limits done by applying pressure or tension is 
known appreciably to change this property, as 
does work which exceeds the elastic limit. Heat 
treatment of material also produces measurable 
effects, as does magnetization of samples. Really 
large changes in e.m.f. from that exhibited by a 
pure metal may be produced by small alloying 
constituents and so on. With these facts in mind 
a form of apparatus which allows the measure- 
ment of the thermoelectric e.m.f. for each incre- 
ment of length of a metal test bar was devised, 
with the expectation that such an arrangement 
would have many useful applications where 
physical changes as a function of position along 
such a bar are to be investigated. The first type 
of apparatus to accomplish this purpose was 
completed a number of years ago and has been 
described in detail elsewhere.! The present ap- 
paratus is an outgrowth of that one, though in 
no sense merely a refined version as there are 
differences in principle as well as detail involved. 

Many obvious uses for an apparatus of this 
kind immediately suggest themselves. For in- 
stance, it could be used to check the uniformity of 
polycrystalline and single crystalline bars and 
thus serve as a production control. A study of 
residual effects resulting from various stress- 
strain conditions applied to a test bar could be 
made both in the elastic and plastic deformation 
regions. Experiments concerning heat-treatment 
and recrystallization effects as a function of some 


' A. Goetz and M. F. Hasler, Phys. Rev. 36, 1752 (1930). 


parameter, such as annealing temperature along 
the test bar, could be readily carried out. Thermo- 
electric effects produced by magnetization of 
various portions of a metallic bar could be in- 
vestigated. In the study of metallic alloys this 
type of analysis is particularly fruitful as the 
function connecting thermoelectric e.m.f. with 
alloy-constituent concentration parallels the 
physical-structive, ‘concentration function. The 
separation of alloying components upon crystal- 
lization also may be studied to advantage by this 
method and made to yield much information 
concerning crystal mechanics.” 


THE PRINCIPLE OF THE THERMOELECTRIC 
ANALYZER 


If a metal bar, placed in a recess possessing two 
regions of constant but different temperature, 
connected by a region having a steep heat grad- 
ient, Fig. 1, has its thermoelectric e.m.f. meas- 


Temperature 


Length 


Fic. 1. Temperature distribution within the thermoelectric 
analyzer. 


2 Much literature exists on these various points; see any 
handbook of physics for a resumé. 
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ured by means of wires fastened to its ends, a 
value will be obtained which corresponds to the 
e.m.f. of that portion of the bar which subtends 
the temperature gradient. If the bar is moved 
through this gradient at such a rate that approx- 
imate thermal equilibrium is maintained, then 
the thermoelectric e.m.f. as a function of position 
along the bar may be recorded for its entire 
length. This constitutes the so-called method of 
thermoelectric metal-crystal analysis. 


METHOD OF PRODUCING THE TEMPERATURE 
DISTRIBUTION 


Since in this thermoelectric analyzer one is 
interested in changes along the length of a metal 
bar as well as the overall thermoelectric e.m.f. of 
bar to copper, say, which is usually measured, a 
very high degree of constancy is required, both 
in the temperature along the regions of constant 
temperature (Fig. 1) and in the difference of 
temperature between these regions. The magni- 
tude of thermal fluctuation allowed was one-tenth 
of a percent of the temperature difference used. 
The method which achieves this result is essen- 
tially a self-regulating one in that differences of 
temperature which occur for any reason occur 
simultaneously in the two regions of different 
temperature and at a relatively slow rate, so that 
along the length of each, constancy of tempera- 
ture is maintained as well as a constant difference 
between the two. Further, in order to maintain a 
steep temperature gradient between the two 
regions a large heat source and sink must be 
maintained. All these various requirements are 
met by using the method of water flows. Two 
identical systems thermally and electrically in- 
sulated were devised, Fig. 2, A and B, and are 
fed from a common source of water. The rate of 
flow through each is identical and since their 
heat capacities are likewise identical, the two 
systems fluctuate in unison with the temperature 
of the common water source. This remains the 
case when one of the flows is heated at a constant 
rate. Since the flows enter the apparatus at dif- 
ferent temperatures, the cold one at the bottom 
of A, the hot one at the top of B, at adjoining 
points, a steep temperature gradient is main- 
tained across this region and, since further heat- 
ing or cooling does not occur, the remainder of the 
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Fic. 2. Schematic diagram to illustrate the principle of 
the thermoelectric analyzer. The water streams keep A and 
B at different temperatures producing the temperature 
distribution along A and B, depicted in Fig. 1. 


regions traversed by the flows are kept at con- 
stant temperatures. Constant pressure is main- 
tained in the water systems by having the water 
fall through a constant height between inlets and 
outlets, while constant heating is effected by an 
electrical heater operated with large capacity 
storage batteries. 


THE THERMOELECTRIC ANALYZER PROPER 


In Fig. 2 chambers A and B are depicted with a 
continuous slot running their full lengths. It is in 
this slot that the metal bar to be analyzed moves. 
The chambers have a cross section as illustrated 
in Fig. 3, which also shows the position of the 
crystal at the center. Since the bar can only make 
contact with two walls its thickness is the only 
critical dimension, as upon it alone depends the 
thermal contact with these walls. In another 
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Fic. 3. A detailed drawing of the chambers A and B through 
which the crystal to be analyzed is moved. 


paper’ a method of growing single crystals pos- 
sessing very exact dimensions was described, as 
was also a means of casting polycrystalline bars 
having this quality. Such test bars are suitable 
for investigation with this type of apparatus. 
The constructional details of the chambers are 
also shown in Fig. 3. The materials used are 
copper and brass, the former being employed at 
all places where the heat conductivity must be a 
maximum, the latter at all other places. Thus 
the cavity in which the crystal moves has walls 
made of rectangular copper pieces a, shaped as 
indicated. The slot width which determines one 
crystal dimension is 0.0530’’ and is accurate to 
0.005” along the whole length. Heavy brass side 


°M. F. Hasler, Rev. Sci. Inst. 4, 656 (1933). 


pieces b, slotted to allow the copper pieces to be 
wedged in tightly, hold the assembly rigid and 
supply necessary heat capacity to the system. 
Brass sides c and end plates of 0.010” copper 
sheet complete the chamber walls. The various 
pieces are soldered together with water tight 
joints. 

From Fig. 2 it is apparent that the temperature 
gradient must be maintained in the crystal by 
applying heat to it at the top of the bottom 
chamber and withdrawing it at the bottom of the 
top chamber. Thus at these points the water flows 
must be most effective. To do this copper-sheet 
baffles direct the incoming water against these 
points. In the body of the chamber, to facilitate 
the maintainence of a constant temperature, 
other copper pieces placed parallel to the flow 
have projections which introduce a stirring ac- 
tion. At the outlets similar arrangements are 
used as at the inlets, primarily to make the 
chambers symmetrical. 

The 4wo chambers are fastened together with 
the aid of cross bars d, which are bridged by 
wooden blocks. Metal tie rods hold these various 
blocks together. In a similar manner the two 
chambers are ‘‘hung up,” so to say, by a similar 
set of blocks mounted toward the top of chamber 
A and supported by the main frame of the ap- 
paratus. Wood is used in this manner of support- 
ing to provide thermal and electrical insulation 
between the two chambers. At the point where 
the chambers press together, a blotting paper 
washer 0.020” in thickness, slotted to allow the 
passage of the crystal, furnishes insulation. 


THE WATER SYSTEMS 


The water systems, besides those portions 
made up of the chambers A and B, consist of 
inlet tanks and tubes, flow pipes and discharge 
tubes. The inlet tanks as shown in the schematic 
diagram (Fig. 2) have three openings each; an in- 
let at the bottom, an outlet a few inches above 
the bottom and an overflow near the top which 
provides a constant “head” for the system. Since 
these tanks are fed from the same water main 
and it is desirable to have the two chambers 1 
and B electrically insulated, the inlet tubes con- 
sist of two eight feet 1/4’ rubber hoses. These 
furnish several hundred thousand ohms resist- 
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ance between the inlet tanks. These tanks, which 
are mounted directly under the ceiling, are con- 
nected to the analyzer chambers A and B by 1,4” 
inside diameter copper tubings which branch at a 
point directly behind the analyzer to make a 
left and right connection to the corresponding 
halves of the chambers. The copper outlet pipes 
branch in a similar manner and discharge into a 
sink through glass tubes six feet in length which 
insulate the chambers A and B electrically. 
Adjustable jets on the ends of these tubes allow 
the water flows through the two systems to be 
made equal. This is an important condition since 
if it is not fulfilled changes in temperature which 
occur in the water source will not effect the two 
systems, which are of equal heat capacity, at the 
same rate. 

The water which enters chamber B must be 
heated to provide the temperature gradient be- 
tween A and B. This is accomplished by an elec- 
tric heater made by winding three sections of 
No. 22 silk covered Advance wire on a four-foot 
glass tubing, which is made part of the connecting 
line from the inlet tank to the analyzer chamber. 
The three sections connected in parallel have a 
resistance of 242. 


THE RECORDING SYSTEM 


To enable a continuous record to be made of 
the thermoelectric e.m.f. as a function of position 
along the crystal bar a photographic recording 
device is used. This is driven by a system of gears, 
worms and a rack, as illustrated in Fig. 4. The 
mechanical ratio is such that one revolution of the 
motor moves the plate 0.00096 in. Since it is es- 
sential that a fixed relationship be maintained 
between the crystal bar to be analyzed and the 
photographic plate, the bar is also moved by the 
same motor. To obtain various photographic re- 
solving powers the crystal may be moved at the 
rates 1, 1/2, 1/4, 1/8 the speed of the photo- 
graphic plate. This is accomplished by a sector 
of a circle which is driven by the same worm 
which moves the plate. The crystal bar may be 
connected by means of a steel tape to various 
radii of the sector, which move at the rates men- 
tioned above. To enable changing from one speed 
to another, the analyzer chambers together with 
their entire water systems can be moved along 
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Fic. 4. Diagram of the recording device by which the 
photographic plate is moved in synchronism with the crystal. 
Ratios of plate to crystal speéd obtainable are 1:1; 2:1; 
4:1;8:1 


on supporting rails so that the slot of the analyzer 
in which the crystal moves is directly below a 
driving radius. 

The actual photographic record consists of the 
deflections of two galvanometers which are 
registered by light beams focussing images of 
vertical slits onto the stationary horizontal slit 
shown before the plate in Fig. 4. 


THe CrystaL BAR MOUNTING 


The crystal is supported and electrical connec- 
tions made to it by means of a crystal holder. 
This is made of a wooden strip just thick enough 
to slide in the analyzer slot. A piece the size of the 
crystal to be used is cut out. Fine glass tubes, 
glued near the edges of the strip, contain the con- 
necting wires which pass through the wood to the 
ends of the opening. When the holder is in the 
analyzer these tubes move in the wider slots at 
the side of the crystal-containing slot as shown in 
Fig. 3. Electrical contacts which also hold the 
crystal bar in place are made by soldering the 
connecting wires to the crystal in little recesses 
at the ends of the holder opening. For low melting 
point crystals, Wood’s metal works quite well as 
a solder. A bent chromel wire electrically heated 
makes a suitable soldering iron for these small 
openings. 

Besides the crystal this holder also carries a 
copper-constantan differential thermocouple to 
measure the temperature difference between the 
two chambers A and B. This is mounted in a 
separate glass tube which runs the length of the 
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Fic. 5. Wiring diagram showing connections to the 
crystal and the differential thermocouple as well as those 
to the various recording and calibrating instruments. 


holder. The constantan portion extends from a 
point adjacent to the center of the crystal bar to 
the lower end of the holder, where it connects with 
the lower crystal contact wire. These various 
connections are made clear in Fig. 5. 


THE ELECTRICAL SYSTEM 


Two galvanometers are used, one to measure 
the crystal bar’s thermoelectric e.m.f., the other 
to measure the differential thermocouple’s e.m_f. 
and hence AT. Since in general the changes in 
e.m.f. along the crystal bar will be small in 
comparison to the overall crystal to copper 
e.m.f., the latter is compensated with a bucking 
circuit so that high galvanometer sensitivities 
may be used in recording the former. This same 
circuit, Fig. 5, is used to calibrate the galvanom- 
eters after each run. This calibration is recorded 
directly on the same photographic plate so that 
all the necessary information for calculating 
absolute values of the thermoelectric power 
AE/AT for any point on the crystal is contained 
in a permanent record. 

The bucking circuit consists of a storage cell, a 
high precision voltmeter, a 1-100,000 ohm re- 
sistance box, and a 1-ohm standard resistance. 

The highest sensitivity obtainable with the 
crystal e.m.f. measuring galvanometer adjusted 
for critical damping was 4.4 microvolts per cm 
at the plate while that used on the other galva- 
nometer was 24 microvolts per cm. The distance 
from the plate to the galvanometers used was 
about 1 meter. 


RESOLUTION OF THE ANALYZER 


The resolution which is obtainable with the 
thermoanalyzer depends upon the temperature 
gradient which can be maintained between 
chamber A and B. This, in turn, depends upon 
the heat conductivity of the bar being analyzed 
and the amount of hot and cold water which 
can be run through the chambers per unit time. 
The first factor, for a bar of given material, de- 
pends only on the thickness of the bar, the 
analyzer slot being built accordingly. The second 
factor, for a given flow resistance determined by 
the dimensions of the apparatus, depends upon 
the non-fluctuating pressure head and power 
supply available. 

Values for these various factors as used in the 
analyzer depicted herein may be of interest. 
The metal bars investigated were bismuth 
single crystals rectangular in shape with the 
critical thickness 0.051’’. This size fits easily into 
the 0.053” slot of the analyzer, the unoccupied 
thickness being filled with light lubricating oil to 
give thermal contact. Bismuth has a low heat 
conductivity so that a steep temperature gradient 
could be maintained along the crystal without 
excessive flows and the correspondingly high 
heating currents necessary to raise the tempera- 
ture of the hot water a few degrees. With 600 
watts on the heater a temperature difference of 
4.75°C exists between the two chambers for a flow 
of 1840 ce per second. This whole temperature 
drop occurs over a 2 mm section of the crystal. 
Thus the thermoelectric e.m.f. for every 2 mm 
piece of crystal length is entirely resolved so to 
speak. Thermoelectric inhomogeneities of much 
less length than this can be partially resolved 
providing they represent large changes in the 
e.m.f. over that of the surrounding regions. 


A SPECIFIC EXAMPLE 


A typical record produced by the thermo- 
electric-analyzer is shown in Fig. 6. The crys- 
tal analyzed therein is a 12 cm bismuth single 
crystal alloyed with 0.04 percent Pb and grown 
at a very slow rate, 0.0003 cm/sec. to illustrate 
how the lead distributes itself within the crystal. 
The wavy line record from left to right gives the 
crystal analysis from front to rear. It is obvious 
from this record that the lead is pushed to the 
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Fic. 6. A typical record obtained with the thermoanal- 
yzer. The wavy line is the thermoelectric analysis of a bis- 
muth crystal, the abscissa representing position along the 
crystal, the ordinate thermoelectric e.m.f. The smooth line 
gives the differential thermocouple e.m.f. 


rear of the crystal by the process of crystallization 
and that increasing lead content raises the ther- 
moelectric e.m.f. for this type of crystal, trigonal 
axis normal to rod. The fine structure on the 
record gives some idea of the length of inhomo- 
geneity which may be partially resolved. The 
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smooth line record is that produced by the difier- 
ential thermocouple and gives the temperature 
distribution as measured by it within the ap- 
paratus. The gradient shown here is not as steep 
as that within the crystal as it is measured’by a 


couple in a glass tube. The temperature differ-. 


ence between chamber A and B is given ac- 
curately, however, by the distance between the 
parallel portions of the record and thus A7 is 
determined. Another method is to measure the 
temperatures of the two emergent water streams 
with a mercury thermometer. 

This record is only one of a large number made 
on alloyed bismuth crystals, which constitutes a 
research problem which will shortly be reported 
in the literature.‘ It serves to illustrate the type 
of thing which may be done to advantage with 
an analyzer of this kind. 

It is a pleasure to thank Dr. Millikan and Dr. 
Goetz, who have made the building of this ap- 
paratus possible. The latter will publish an article 
shortly on specific uses for an apparatus of the 
kind described herein. 


*To be published shortly in the Physical Review, if 
possible. 
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The Shot Effect and Electrical Breakdown in Insulators! 


F. E. HAwortH AND R. M. Bozortn, Bell Telephone Laboratories 
(Received October 26, 1933) 


Measurements are made of the fluctuations in currents 
passing through insulating materials subjected to high, 
constant electric fields. Applying the theory of the shot 
effect to these fluctuations we are able to calculate the 
average number of elementary charges (electrons or ions) 
which discharge onto the electrodes as a unit. In glass, 
this number increases exponentially with the field from 
about one at low fields to about 10° at fields just below 
breakdown. Considered with respect to various theories 


offered in explanation of electrical breakdown, our results 
indicate that the shot effect is caused by ‘“‘avalanches”’ of 
electrons discharging through paths in the material over 
which breakdown occurs at higher fields, the avalanches 
increasing in size with the field until the material is 
ruptured. By making use of the Johnson effect (thermal 
agitation of electricity) a new and simple method has been 
worked out for calibrating the amplifier for shot effect 
measurements. 


INTRODUCTION 


INCE Schottky? first discovered and inter- 

preted the shot effect in vacuum tube circuits, 
the theory has been applied to a variety of prob- 
lems. The fact emerges that studies of the elec- 
trical fluctuations or “‘noise’’ accompanying 
many physical phenomena, amplified as may be 
necessary, are valuable in getting a better insight 
into the physical processes involved. We have 
recently measured the fluctuations in currents 
through insulating materials subjected to high 
electric fields, and as a result have been led to 
certain conclusions regarding the nature of elec- 
trical breakdown in solids. 


METHOD OF MEASUREMENT 


The apparatus is shown diagrammatically in 
Fig. 1 and was similar to that usually used for 
measurements of the shot effect. The insulator to 
be tested was at D between the plates of a small 
condenser. The battery B forced through it a 
current which also passed through R, the input 
resistance of the amplifier, so that current fluctua- 
tions in the insulator impressed voltage fluctua- 


tions upon the grid of the first tube. The ampli- 


' Presented in part before the Committee on Electrical 
Insulation of the Division of Engineering and Industrial 
Research of the National Research Council, Boston, 
November 13, 1931. Also Phys. Rev. 39, 845 (1932). 

?W. Schottky, Ann. d. Physik 57, 541 (1918); 68, 157 
(1922). 
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fied fluctuations were impressed at the output 
upon a vacuum thermocouple connected to a 
galvanometer, the deflections of the latter thus 
measuring the fluctuations at the input of the 
amplifier. By a switching arrangement the same 
galvanometer was used to measure the magnitude 
of the conduction current in the insulator. In a 
telephone receiver substituted for the thermo- 
couple, the fluctuations could be heard as a 
sharp crackling noise. 


Ao 


Fic. 1. Noise measuring circuit. D is the specimen, B a 
battery variable from 0 to 1500 volts, R a resistance of 
0.1 to 0.5 megohm, Th a thermocouple and G a galvanome- 
ter. 


The amplifier had a voltage amplification fac- 
tor of about 10°, and since the input circuit was of 
very high impedance, several precautions were 
necessary to eliminate all sources of noise outside 
the specimen. The battery B was made up of dry 
cells to a maximum potential of 1500 volts and 
had to be very carefully insulated to eliminate 
small leakage currents. The amplifier and the 
whole input circuit were carefully shielded. 
Furthermore, a grounded guard ring or other 
suitable means was provided on the specimens to 
eliminate the effects of surface conduction. 

The insulator used in most of the experiments 
was G702EJ Pyrex glass, although measurements 
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were also made on lead glass and mica. The glass 
specimens were in the form of bubbles blown from 
small tubes. These little flasks were filled with 
mercury and immersed in a small beaker of 
mercury, the mercury surfaces inside and outside 
being at approximately the same level. Since the 
edges of the mercury electrodes were separated 
by glass much thicker than the bubble, the so- 
called edge effect was eliminated.’ The mica 
specimens were in the form of flat sheets, upon 
which silver was chemically deposited to form 
electrodes, and these were placed between flat 
plates with carefully rounded edges. 

Since the fluctuations due to the thermal agita- 
tion’ of electrons in the resistance R could not be 
eliminated, they were allowed for in the measure- 
ments. Their presence is actually advantageous 
since it permits us to eliminate, as will be shown 
below, factors which are difficult to measure. 
It was found that a wire-wound resistance or one 
made of a mixture of xylene and alcohol was 
suitable for the resistance R, but that an ordinary 
grid leak could not be used on account of addi- 
tional voltagé fluctuations produced by the pas- 
sage of current through it. 

Several final tests were made to prove that the 
fluctuations originated in the specimen. The in- 
sulator under test was removed and all parts of 
the circuit so thoroughly insulated that no leak- 
age current could be detected either directly or 
by the fluctuations such a current would produce. 
In a second test a resistance of the same order of 
magnitude as the specimen, but known to give no 
additional fluctuations upon the passage of cur- 
rent through it, was substituted for the specimen 
whereupon the fluctuations ceased. They also 
ceased when one-half the battery was reversed, 
the potential applied to the insulator being thus 
reduced to zero while all the other elements of 
the circuit remained the same. Under these condi- 
tions if noise had originated in fluctuations of 
battery potentials, or in electrical pick-up due to 
insufficient shielding, it would still have been 
present. 

In order to show that the noise arising in the in- 
sulator is not due to electromotive forces arising 


3L. Inge and A. Walther, Archiv f. Electrotechnik 19, 
257 (1928). 
‘J. B. Johnson, Phys. Rev. 32, 97 (1928). 
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from local changes in chemical composition, we 
used a Pyrex glass bulb with a 90 percent sodium 
amalgam as electrode material in place of mer- 
cury. The flask was baked out in a high vacuum 
and the amalgam introduced without coming in 
contact with air. With this arrangement any 
metallic ions plated out of the glass are replaced 
by ions of the same sort. Also, there is no question 
here of a discharge taking place in a thin film of 
gas between the electrode and the glass. We found 


no difference in the fluctuations due to this. 


change in the electrodes. 


THEORY 


From the power developed in the measuring 
circuit, we can calculate a quantity which de- 
pends upon the sizes of the groups of charges 
which arrive as a unit at one of the electrodes. 
The mathematical theory of the shot effect has 
been worked out in detail, and we follow the 
presentation of Fry® with certain extensions for 
our special circuit. He derives: 

S=) 
where S is the power dissipated in the measuring 
circuit (in our case in the thermocouple) as a 
result of the fluctuations, 7 is the mean frequency 
of the electrical impulses resulting from the ar- 
rival of groups of charges at the electrodes, and 
w is the mean energy dissipated in the measuring 
circuit for one impulse. For W, he derives the 


relation: 
0 


where ” is the number of charges in the group, 
e is the electronic charge, and E(t) and i(t) are 
the time-functions of the voltage and current in 
the measuring circuit due to unit charge placed 
on the electrode at time /=0. 

But ip=v(ne), where iy is the mean current 
through the specimen as measured by a direct- 
current instrument, giving 


(me)*> 7 
] D(t)a(t)dt. 


5’ T. C. Fry, J. Frank. Inst. 199, 203 (1925), 
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Simplified by the Carson-MacColl theorem® 
and applied to our special circuit this reduces to: 


wr 


Here ip is the current in the thermocouple, R is 
the input resistance, 7 is the transfer admittance 
of the amplifier, w is 27 times the frequency of the 
electrical oscillators and C is the total capacity 
across the input of the amplifier, including the 
tube capacity, that of the condenser formed by 
the electrodes and the specimen and the capaci- 
ties of the connecting wires. 

Instead of evaluating this integral by calibrat- 
ing the amplifier and integrating the resultant 
curve graphically, we observe that for this circuit 
Johnson! obtained the relation: 


where 7 is the current produced in the thermo- 
couple by the thermal agitation of electrons in 
the resistance R, k is Boltzmann’s constant, T 
the absolute temperature and C has the same 
meaning as in the previous equation. Dividing 
ip? by ir*, we have the relation: 


= 
which is independent of C and of the character- 


istics of the amplifier. In our circuit whenever ip 
is present, 77 is also, so that we actually measure 


ip?+t7*; but by reducing the voltage on the insu- . 


lator to zero, we measure 77” and can therefore 
obtain ip? by subtraction. 

Before using the relation obtained above we 
tested it by measuring the shot effect in a two 
element vacuum tube which was substituted for 
the insulator. For such a tube 7 is one and there- 
fore n2/n is one, and it should be possible to use 
the circuit for an accurate measurement of e/k. 
By using a Western Electric No. 246A vacuum 
tube as the first amplifier tube and a Davohm 
wire-wound resistance, it was found that if R was 
kept below half a megohm the expected value of 
e/k was precisely that obtained by measurement. 
In this experiment 7) was much larger than in the 
case of insulators and it was necessary to keep it 
so low that the potential drop through R had a 
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negligible effect upon the grid of the first ampli- 
fier tube. We also found that it was necessary to 
put a high pass filter in the output circuit with a 
cut-off at about eighty cycles, so that the “‘flicker 
effect’’® and all remnants of sixty cycle pick-up 
were effectively eliminated. 

This method suggests itself as a possible one 
for determining the ratio e/k with precision. But 
this ratio is already known with much greater 
precision than the separate values of e and k and 
is in fact equal to F/R, the Faraday divided by 
the gas constant. With this method it is not to 
be expected that the errors could be kept below 
the present probable error in e k which is about 
0.01 percent. 

As a further check upon our circuit we cali- 
brated the amplifier by means of an oscillator, 
determining yy, as a function of the frequency. 
This was plotted and integrated graphically and 
by use of the result, e was measured from the shot 
effect in the vacuum tube and k from the resist- 
ance noise in R, these values being in substantial 
agreement with the known values. 


RESULTS 


Fig. 2 shows a plot made from typical noise 
data for Pyrex glass. As may be seen, these data 
scatter badly. However, the fluctuations always 
increased with increasing potential, and by mak- 
ing many different runs on a given specimen, a 
good average value for the curve was obtained. 
Fig. 3 shows the conduction current 7p as a func- 
tion of the applied voltage for the sample of 
Pyrex used for the noise measurements plotted in 
Fig. 2 and in arbitrary units the conductivity of 
this specimen as dependent upon voltage. 

The quantity ?/n is a kind of average value of 
the number x of unit charges in a single group, 
and attains values as high as 10° in glass. These 
experimental values may be compared with the 
values of predicted by certain theories to be 
discussed in the next section. 


SIGNIFICANCE FOR BREAKDOWN 


The experimental results may be considered in 
relation to various current theories of electrical 
breakdown. In 1927 it was suggested separately 


6 J. B. Johnson, Phys. Rev. 26, 71 (1925). 
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Fic. 2. Group-size as dependent upon voltage. Here n 
represents the number of charges equal in magnitude to 
the charge of one electron which discharge onto the 
electrodes together. The quantity determined by experi- 


ments on the shot effect is the ratio of the means, n?/n. 


by Joffé,? Rogowski® and Smekal® that in rock- 
salt, as an example of an insulator, breakdown 
came about as the result of an excessive positive 
ion current. According to their views, the mobile 
Na* ions were speeded up by the action of the 
field until they acquired sufficient energy to free 
more ions by impact. In Joffé’s theory the ions 
were supposed to be loosened from their normal 
positions in the crystal lattice; in the theories of 
the others, the new ions freed were those loosely 
bound in imperfect parts of the crystal, the im- 
perfect region being enlarged as the ion-current 
and field-strength increased. 

More recently Rogowski,'® von Hippel! and 


7A. Joffé, T. Kurchatoff and K. Sinelnikoff, J. Math. 
and Physics (M.I.T.) 6, 133 (1927). A. Joffé, The Physics 
of Crystals, McGraw-Hill, New York, 1928. 

8 \W. Rogowski, Archiv f. Elektrotechnik 18, 123 (1927). 

® A. Smekal, Archiv f. Elektrotechnik 18, 525 (1927). 

10 W. Rogowski, Archiv f. Elektrotechnik 23, 569 (1930). 

A. vy. Hippel, Zeits. f. Physik 67, 707 (1931). 
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Fic. 3. Conductivity and current as dependent upon 
ee for the same glass bulb as used for the data of 
ig. 2. 


others” have proposed that the breakdown is 
produced by an “‘avalanche”’ of electrons due to 
acceleration by the field and multiplication of the 
number by ionization by impact. The three chief 
experimental facts supporting the electron theory 
as contrasted with the ion theory are: (1) the 
short time (10-7 sec.) of breakdown" considered 
with reference to the low mobility of the ions as 
compared with the electrons; (2) the lack of 
dependence of breakdown potential upon the 
temperature’ and therefore on the ion current 
(the resistivity” of rocksalt increases by many 
powers of ten as the temperature decreases from 
that of the room to that of liquid air, whereas the 
breakdown voltage does not then increase'®); (3) 
the Lichtenberg patterns obtained in the solid 
by von Hippel, indicating that cumulative ioniza- 
tion is caused by negative carriers. 

Our results are readily interpreted in terms of 
the electronic theory. According to this theory 
the electrons make a negligible contribution to 
the ordinary conductivity’? even at relatively 
high fields, but when the field is sufficiently high 
and the ionizations by impact sufficiently great in 
number a channel is formed from one electrode to 


2, Rosenkewitsch and C. Sinelnikow, Zeits. f. Physik 
73, 118 (1931); W. Werner, Ann. d. Physik 86, 95 (1928). 

8 WW. Rogowski, E. Flegler and R. Tamm, Archiv f. 
Elektrotechnik 18, +79 (1927). 

44 \WV. Werner, reference 12; A. v. Hippel, reference 11; 
H. Rochow, Archiv f. Electrotechnik 14, 361 (1925). 

' A. Joffé, reference 7 (book), p. 129. 

© W. Werner, reference 14. 

17H, Schiller, Ann. d. Physik 83, 137 (1927). 
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the other in much the same fashion as in a dense 
gas. We should expect the average number of 
electrons in such an avalanche or group of 
charges to increase exponentially with the field. 
No very definite statement can be made about the 
expected number of electrons per avalanche at 
breakdown; but Rogowski'* in considering the 
possibility of observing a single group of charges 
on a cathode-ray oscillograph has estimated that 
in air at atmospheric pressure with electrodes 1 
cm apart the number of unit charges per group 
may be of the order of 10%. We may also estimate 
that at breakdown the number in a group will be 
larger by several orders of magnitude than the 
number of atoms in a single row reaching from one 
electrode to the other, for at breakdown we ex- 
pect a channel to be formed within which the lat- 
tice is disrupted. In the thin glass used in our ex- 
periments the number of atoms in a single row is 
about 10°. 

The expected law of increase of group size with 
voltage is confirmed by our experiment, and the 
above rough estimates of size are at least not in- 
consistent with our results. On the other hand, 
the group-size predicted by Joffé’s theory is hope- 
lessly small compared to the experimental values. 
The number per group predicted on his theory is 
equal to the ratio of the conductivity at the given 
field to the limiting conductivity at very low 
fields, and in the case of our glass specimen is less 
than 10 as shown in Fig. 3." 

To express it in another way: if the increase of 


'8 W. Rogowski, Archiv f. Electrotechnik 25, 596 (1931). 

'? Recent experiments of G. Just, Zeits. f. Physik 82, 
119 (1933), even throw considerable doubt on the validity 
of Joffé’s determinations of the breakdown of very thin 
films and indicate that the high values he obtained were 
due to adsorbed films. It was on the basis of his own 
results on thin films that Joffé was led to his theory of 
ionic breakdown. 
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conductivity in high fields is due to an increase in 
the number of mobile positive ions produced by 
ionization by impact, the “‘noise’’ we should 
expect from this source alone is orders of magni- 
tude less than that which we actually observe at 
the highest fields. If the greater part of the noise 
is caused by the presence of relatively few groups 
of electrons, the average population of these large 
groups cannot be determined unless we know the 
amount of the current carried by the electrons 
alone. Tests of Faraday’s law of electrolytic con- 
duction” in similar materials indicate that in 
high fields the positive ions account for at least 
95 to 99 percent of all of the current, so that we 
may set the upper limit of the fraction carried by 
electrons as about 10-?. The average number 7, 
for the electron groups corresponding to this frac- 
tion is then 10? times the previously calculated 
value of 2?/n, or about 10% for our material at 
fields just below breakdown, as obtained from the 
final equation above by substituting for ip the 
strength of the electron current. The passage of 
this number of electrons in a single avalanche 
recalls the experiments of Inge and Walther*! and 
of von Hippel'' who observed microscopic chan- 
nels in rocksalt and glass after a field-strength 
ordinarily sufficient for breakdown was applied 
for a very short time. It was only after the field 
was applied for a longer time that the channels 
were large enough to cause breakdown. 

We may conclude that at the present time the 
experimental data, including those reported here, 
point to a mechanism of breakdown by electron 
avalanches, and that the electrical fluctuations 
which we have observed are due to smaller 
avalanches insufficient for breakdown. 


20H. Schiller, reference 17; A. Joffé, Phys. Zeits. d. 
Sow. 1, 176 (1932). 


21 L. Inge and A. Walther, Archiv f. Electrotechnik 24, 
259 (1930). 
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This paper is an analysis of the unsteady flow of a com- 
pressible fluid flowing radially to a well in a sand forma- 
tion. The phenomenon of unsteady flow occurs as a result 
of fluid expansion. When the pressure in the formation is 
lowered, the fluids therein expand and the increase in vol- 
ume imparts motion to the fluid which flows towards the 
region of lowest pressure in the formation. This process is 
continuous in the reservoir and extends further away from 
the well with increased production. By the “equation of 
continuity,’’ the solution for two specific cases of unsteady 
flow are derived. The first of these is that in which the fluid 
from a sand reservoir of limited size flows to a well in 
which the pressure at the level of the producing sand always 
remains constant. The variation of the pressure gradient, 
the rate of production, and the cumulative production with 


respect to time are given in Eqs. (11), (24), (29) and (30). 


The second case is that in which the flow of fluid to a well 
is such that the rate of production at the well is always con- 
stant. This case is derived by the assumption that the 
cylindrical body of sand, which is influenced by the well, is 
subject to a steady depletion of fluid, and in order that the 
rate remain constant at the well, fluid must flow into the 
reservoir from some extraneous source in increasing 
amounts. Eventually, however, the rate of flow of fluid 
into the reservoir becomes equal to the rate at which fluid 
is withdrawn from the well and steady flow is established 
in the sand. The equations for the pressure gradient and the 
decline of pressure at the well with respect to time are (39), 
(52) and (56). 


INTRODUCTION 


LUIDS encountered in oil reservoirs are 
compressible, and because of their expansion 
during the course of production, the flow of these 
fluids to wells is unsteady. When the well is 
opened, the pressure at the well bore is lowered, 
causing: the fluid near the well to expand and 
flow toward the well. The loss of fluid from the 
reservoir lowers the pressure in the more distant 
portion of the reservoir, permitting the fluid in 
those portions to expand and flow to the well. 
Thus, the area around the well affected by drain- 
age is continuously enlarging, and the pressure 
and velocity of the fluid is changing continuously 
at all points in the reservoir. This condition is 
characteristic of the unsteady state of flow. 
Under certain circumstances, the equations for 
steady flow are useful approximations. This is 
true because the rate of change of the pressure 
gradients in the later periods of production is 
relatively small. The analysis of steady flow of 
fluids through sand formations has been de- 
veloped by different investigators.': 


1Schlichter, U. S. G. S., Nineteenth Annual Report, 
1897-8. 
? Morris Muskat, Physics 2, 329 (1932). 
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EQUATION OF CONTINUITY 


In order to develop the analysis for the un- 
steady flow of fluids through a sand, certain as- 
sumptions must be made. These are as follows: 


(1) The permeability, porosity, and thickness of a sand 
are uniform throughout that portion of the formation 
affected by a well. 

(2) The well bore is cylindrical and penetrates the entire 
thickness of sand. 

(3) The mixture of oil and gas behaves as a homogeneous 
fluid, and the velocity of flow at any point in the sand is 
directly proportional to the pressure drop per unit of length 
along the path of flow. 

(4) The density of the fluid within the sand is expressed 
by the equation 
e( P— Pw) (1) 
where py is the density of the fluid at the well, P, the pres- 
sure at the well face, and c the compressibility of the fluid. 


P= Pw€ 


The exponential expression for density is an 
assumption which has been adopted for con- 
venience. It is applicable for a slightly compres- 
sible fluid‘ since it reduces to 


p= Pwll +c(P—P.) ], 


3H. D. Wilde, Jr., and T. V. Moore, Oil Weekly, Decem- 
ber 5, 1932. 

‘Lamb, Hydrodynamics, Cambridge University Press, 
p. 6, 1916. 
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solubility and gas laws; (2) calculated from equation 
p=pye?—Pw), 


when c is small. Furthermore, for a fluid of rela- 
tively large compressibility, it gives a fairly close 
approximation over a wide pressure range as 
shown in Fig. 1.5 

The equation of continuity, which is _pri- 
marily a fluid balance over an infinitesimal ele- 
ment of sand, is expressed by the equation 


O(pu) A(pv) 
= 
Ox oy 


(2) 


where u and v are the component velocities of the 
fluid which flows parallel to the plane YOY and 
toward the origin, Fig. 2. The constant, ¢/¢’, is 
the ratio of the volume porosity of sand to its 
surface porosity and @ is time. 
The component velocities of the fluid may be 
expressed as 
K oP 
u=——— and 
Ox 


K oP 
(3) 
¢’ dy 
where K is the permeability of the sand. 
Therefore, 


pu= — p(K/$')(dP/dx), 


but 
0p/dx AP /dx=cpdP, dx; 
then 
K dp dp 
and w= -——. 
Ox dy 


The substitution of these values in Eq. (2) 
gives 
(0° a? =dp/ 00, (4) 


where a? = K/c¢. 


Because this equation is similar to that for the 


° Moore, Schilthuis and Hurst, Oil Weekly, 69, May 21, 


1933, 


IK 


unsteady flow of heat through a homogeneous 
solid, a? may be referred to as the ‘hydraulic 
diffusivity,’’ as comparable to the ‘‘thermal dif- 
fusivity” in the heat equation.°® 

For the radial flow of fluid to the well, this 
equation is transformed so that 


00 (5) 


and the solution of this equation is 


0(Brr) +B, Vo(Bpr) (6) 
n=l 
In this expansion Jo(8,7) and Yo(8,7) are 
Bessel functions of zero order. The second kind of 
Bessel, Yo(8,7), is the Weber-Schlafi? function. 
Furthermore, A,, B,, and 8, in Eq. (6) are con- 
stants to be determined from the conditions im- 
posed on the problem. The possibility of this 
expansion’ and that the series can be integrated 
term by term is assumed for the present. 
Associated with these Bessel functions are the 
recurrence formulae® which are applicable in this 
analysis. 


(7) 
d\2-"J,(2)} —2"Jn4i(2), (8) 
d{\z"Y,(z)} /dz=2"Y,_1(2), (9) 


6 Carslaw, The Conduction of Heat, Macmillan, p. 8, 
1921. 


7 Numerical tables for both kinds of Bessel functions are 
given by Watson, Theory of Bessel Functions, Camb. 
University Press, 1922. 

8 Hobson, London Proc. Math. Soc., [2] 7 (1909). 

’ Watson, reference 7, 2.12 and 3.56. 
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Case I. CONSTANT WELL PRESSURE 

lhe limits for the case when the well pressure remains constant and the depletion of fluid takes 
place in a finite cylindrical sand reservoir, concentric with the well, are as follows: r 
r=a, P=P, 


O<A<x 
r=b, lim. (@P/ar) 


and 6=0, P= Pr. From Eq. (1), these conditions may be interpreted in terms of density, and a plot 
illustrating these limits is shown in Fig. 3. It is observed that the field is limited in extent by the as- 
sumption that the pressure slope is zero at the radius 0, and therefore no fluid can flow across this 
boundary. 
Let 
p=(pr—pw)At+ pw, (11) 


where p,, and pr are constants. Then by Eqs. (5) and (6) 
A= [A nJo(Bnt) +Bn Yo(Bnr) (12) 


As a simplification of this problem, let 


Vo(Bnr) = Yo(Bna@)Jo(Bnr) — Jo(Bna) Yo(Bnr), (13) 
then 


Mn (14) 


where A,,=./,, Yo(8,a) and B,= — \/,J(8,2). Therefore, 
lim. A=¥ My = 0 
ra n=l 


and then by Eq. (11) p=p,, when r=a. 
The differentiation of the series expansion, Eq. (14), with respect to r is the equation 


= Ba MnV (Bur 


n=l 


10 
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where 
V => } —Jo(B,a) J (15) 
by the recurrence formulae Eqs. (8) and (10). If 
OA 
lim. —= = 0 
rb Or n=1 


the limit at the outer boundary 6 is satisfied if we assume that each term in the series expansion is 
equal to zero. Then, 1\(6,5) is equal to zero or 


Jo(Bra) Y,(8,a) Ji(B,b) Y,(8,)) (16) 


and Eq. (16) establishes the many values of the root 8, in the expansion from n=1 to n= <=. 
Finally, in order to establish the limit when time is zero it is necessary that 


lim. A= MV + MeV ol Bor) + MyVo(Bsr) (17) 


=1 
whena<r=b. 
The coefficients 1/,, \/2, etc., may be determined by the following method: 
Let I’(rr) and Io(wr) be any two Bessel functions whose roots are not equal. Then, 
rd? Vo(rr)} /dr’+rd {| Vo(rr)} /dr+72r?Vo(rr) =0 
(1/r)(d/dr) (rd Vo(rr)} /dr]+7?Vo(rr) =0. 


(1/r)(d/dr) [rd { Vo(wr)} /dr]+2Vo(wr) =0. 


or 


Likewise, 


By combining these two equations, we have 


df d\Vo(rr)} Vo(wr)} 
(w? — rr) Vo(wr) = ~|- 
dr dr d dr 


The integration of both sides of this equation with respect to 7 from the limits a to b gives 


d}Vo(rr)} d{ Vo(wr)} 


dr dr 


a 


Since (7b), Vi(wb), Vo(ra), and Vo(wa) are equal to zero by the conditions already established then, 


b 
f rVo(rr) Vo(wr)dr =0 (18) 


a 


when 
If, however, the two roots are equal, this integral may be established as follows: 


(1/r)(d/dr) [rd { Vo(rr)} /dr]+7?2Vo(rr) =0. 
This equation multiplied by Vo(rr)} /dr gives 
(d/dr)(rd Vo(rr)} /dr P+ (d | Vi?(rr)} /dr) =0. 
The integration of this equation with respect to r from the limits a - bis 


6 


b 
rV = (1/27?) (rd {| Vol rr) } /dr)?+(rr Vo(rr))? ] 


a 


a 


=43{[bVo(rb) P—[aVi(ra) 


‘ ta 
yt 
| | 
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since V,(rb) and Vo(ra) are equal to zero. However, 


Vi(ra) = Yo(ra)Ji(7a) — Jo(ra) Y;(ra) 
is comparable to 
Yn(2)In41(2) — Jn (2) Ynii(2) =2/nz (19) 


given by Watson.'® Therefore, V\(ra) =2/a7a and 


f rV = (1/27? x?) { [arb Vo( rb) (20) 


a 


Another definite integral 
b 


f rVo(rr)dr = (21) 


a 


may be derived by the recurrence formulae and Eq. (19). 
If then Eq. (17) is multiplied through by rVo(8ar) where 8, is a root appearing in the series, and 
each term is integrated with respect to r from the limits a to 3, 


b b 
by virtue of Eq. (18). : 


Furthermore, by the substitution of Eqs. (20) and (21) in this formula 


(22) 
and 
78 
A=4r>0 (23) 


n=1 14 [ V0(B,0) F} 


If, however, we adopt the well radius, a, as the unit of length in Eq. (23), and let 8, =8,'/a, r’ =r. a, 
and b’=b/a this equation becomes 


where Vo(Bn'b’) and Jo(Bn’)/ Yo(Bn’) =Ji(Bn’b’)/ Yi(Bn'b’). (25) 


Production rate 


Production rate is equal to the product of the velocity of fluid entering the well and the area of the 
surface of the well bore, or 


dw = 27aH?¢' pa, (26) 
where H/ is the thickness of sand, and yu, the radial velocity of the fluid at the well. Further, 
K K Op 0A 


— lim. lim. — and lim. — =(pr— pw) lim. —. 
Or  Cpy Or r>a OF ra OF 


If we assume that each term in the series expansion Eq. (23) may be differentiated with respect 
to r in the interval a to b, then 


10 Watson reference 7, p. 77. 
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or n=1 ra or n=l 


since 1:(3,a) =2/78,a. Furthermore by substituting 6, =6,’/a then 


dA 2 
lim. — = —— 
Or Td n=1 
and 
= —[4HK(pr— pw)/Cpw] My'e (27) 


n=1 


The rate of fluid produced at the well can be determined as the volume it would occupy at the initial 


) formation pressure instead of well pressure, since 
dwPw = WRPR- (28) 
Therefore, 
= —4 (29) 
KIT pr- pw) n=l 
or 


The reciprocal of this argument versus time is represented as a family of curves in Fig. 4. 
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Cumulative production 


Since cumulative production is the integration of the rate with respect to time, it is expressed as 


6 
Or -f 
0 


Therefore, the integration of Eq. (27) with respect to time and the substitution of (28) gives the 
formula 


Qrprorc co M,/(1 — a?) 


=-4> (30) 
pw) n=1 Bn” 


or 
= F,(b/a, 026/02). 


The cumulative production may be determined graphically when the pressure gradient is estab- 
lished in the sand. This is given by the equation 


Qr b 
=— r(1—A)dr. (31) 
] ad, 


Furthermore, by means of this equation, the total cumulative production of fluid from a reservoir 
over an infinite time can be calculated. 


Fig. 5 
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Since lim. (1 —A) =1 therefore, 


642 


Or 2r 
-— f rdr = —1), 
"(PR Pw] a? dy 


which is the asymptote that the corresponding b’ curve in Fig. 5 approaches as a limit. 


Case II. CONSTANT PRODUCTION RATE 


In order to maintain constant rate of flow at a well, it is necessary that the slope of the pressure 
gradient at the well face remain constant at all times. 

This case is derived by the assumption that the unsteady flow of fluids in a reservoir to a well 
depletes the sand of fluid to such an extent that in order to maintain the rate constant at the well, 
fluids must flow across the outer boundary into the reservoir in increasing amounts. Eventually, 
however, this rate of flow of fluid into the reservoir is the same as the rate of production from the 
well and steady flow is established in the sand. 

Let 

(32) 


where p,, is the density of the fluid at the well in infinite time, and it is assumed that 
A=r, 0)+n(r). (33) 
40, 


Since 


therefore 


(34) 
and 
(35) 


The solution for (34) is as before, 
E= + Bn Yo(rnr) Jer"? (36) 
n=l 
and for (35), it may be shown to be, 
n=Clinr+D. (37) 
This equation expresses the conditions of steady flow since the limit of £ is zero when time approaches 
infinity. 
The limits for this problem are as follows: 


lim. =constant 
or 0=0=« and A=1, a=r=b. 
lim. A=1 | 


The gradients for this case are shown in Fig. 6. 
Steady state 


If we substitute the limits for the steady state in Eq. (32), namely, 


r=da, p=p,, and r=b, p=pr, 
then 


pPr=(pr—p,.)(C In b+D)+p,, and Poo = (pr—p..)(C In a+D)+ p,. 
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Therefore, 
n=In (r/a)/In (b/a). (38) 
The rate of production at infinite time is the equation 
2rallK Op 2rallIK 
Veo = lim. — or qr=- lim. — 
Peo ra Or roa Or 


if the fluid be measured at the original reservoir pressure. However, 


lim. (0/87) = (px — Po) lim. (dn/dr) = (px — In (b/a) 


and 
Po = PR— [arpre In (b/a) |/(27HK). 


Therefore, if p,, is substituted in Eq. (32), 


p= pr—L(1—A)qrere In (b/a) ]/(27HK). (39) 


Unsteady state 


The necessary limits for the unsteady state solution of this problem are as follows: 


Since 
(33) 
and 
lim. —=lim. —+lim. —, then lim. —- =lim. —= constant 
roa or rea or rea or roo or ra or 
if 
ag 
lim. —=0, (40) 
Ta or 
which is the requirement that the pressure slope at all times is constant at the well. 
Furthermore, since 
lim. A=1, and lim. n=1, 
therefore 
lim. (41) 


rb 


| 
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Finally, since A=1 when time is zero, then 


lim. £=1—In (r/a)/In (b/a) = —In (7/6) In (b/a). (42) 


If we let 
= Vo(And) Jo(Anr) Jo(Anb) Yo( An?) 
in Eq. (36), then 


£= (43) 


where A, = N,,Yo(A,5) and B, = — N,,Jo(A,.b). Therefore the limit of £ as ry approaches 6 is zero and the 
condition of (41) is satisfied. 
Furthermore, 


Uy "8, where Uy(rnr) = “o( And) (An?) Jo(And) 


n=1 


by recurrence formulae (8) and (10). 


Then, if 
dg 
lim. —= (44) 
ra or n=1 


and we assume that each term in the series expansion is equal to zero, the condition of (40) is ful- 
filled if U,(A,a) =0. Or, 
Ji Qna)/ ViQana) = Jo(Anb)/ Yo(A,0). (45) 


This equation establishes the values of \, from n=1 ton=<. 
For time equal to zero 
—In (r/b)/In (b/a) = Uo(Ayr) +Ne Uo(ror) +N; Uo(A3r) (46) 


The values of N,, Ne, etc., may be derived by a method similar to the one used in the preceding case. 
It can be shown that 


rU Uo(wr)dr=0 (47) 


a 


when and 
b 


= (1/2772?) (ra) (48) 


a 


when the roots are equal. - 
Another definite integral, which is applicable in this analysis, is 


b 
In (r/b)Uo(rr)rdr = — (49) 


This integration is derived by the recurrence formulae and the conditions imposed on the problem. 
If Eq. (46) is multiplied through by rUo(Anr), where \, is a root appearing in the series, and each 
term integrated with respect to r from the limits a to b, then 


1 Si In (7/b)U o(dnr)rdr 
In (b/a) Seer U (rnr)dr 


by virtue of Eq. (47). By substituting (48) and (49) in this formula 
Ny = 22?U (Ana) /In (b/a) o(Ana) (50) 


co = 
| 
—_ 
) 
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The final form is the expression - 


t= 
In (b/a) »=1 


and 


A=. In (r/a)/In (b/a). (52) 


n=l 


The variation of well pressure versus time may be determined by Eqs. (39) and (52). Therefore, 


Pw=pr— lim. (1—A). (53) 
However, 
1= (r/a)/In (b/a) (54) 
n=1 


for a=r=b, when @=0. So by combining Eqs. (52) and (54) for corresponding roots, we have 


lim. (1 —A) = Nn Uo( Ana) (1 — | 


ra n=1 


If as in the preceding case, the well radius, a, is taken as a unit length, then 

lim. (1—A) = yk (55) 
ra In (b’) n=1 [ An’ ol dn’) 


= Jo(An’b’) / Vo( dnb’). 
Therefore, by substituting Eq. (55) in (53) and simplifying 


=r> (56) 
Gre n=1 hn’) P} 
and P,, is now the variable well pressure. 
The decline curve for well pressure versus time is shown for a particular sand reservoir in Fig. 7. 
The author is indebted to W. K. Lewis, of the Department of Chemical Engineering, Massachusetts 
Institute of Technology, and H. D. Wilde, Jr., in charge of production research at the Humble Oil & 
Refining Company, for the valuable suggestions which they have contributed during the course of this 
work. Furthermore, the author wishes to acknowledge the cooperation of T. V. Moore and R. J. 
Schilthuis whose kind assistance has made possible the application of this theory® to the oil industry. 
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A Magneto-Optic Method of Determining the Vitamin Content of Various Substances 


G. M. Wissink, Physics Laboratory, Iowa State College, Ames, Iowa 
(Received November 14, 1933) 


In this investigation, the magneto-optic apparatus was 
used, but it was modified to some extent from the original 
apparatus used by Allison. The modifications included a 
mercury vapor arc and a monochromator to furnish a 
steady source of monochromatic light, and a photoelectric 
cell and amplifier circuit to read the minima. The apparatus 


gave a distinct minimum, at 32.9 on the electrical path 
scale, for all substances tested which contained vitamin A 
but did not show this minimum for similar substances in 
which this vitamin was lacking. Pure carotene did not give 
this minimum. Pure carotene after irradiation with ultra- 
violet light, however, did give it. 


INTRODUCTION 


VER since the discovery of vitamins, much 
research has been done in an effort to isolate 
the vitamin from the food containing it and to 
find the vitamin content of various foods. In 
order to determine the vitamin content of foods, 
the biologists have used the feeding test method 
but this method has been rather slow and hard 
to control accurately. During the last few years, 
attempts have been made to find a more rapid 
method than the biological for finding this quan- 
tity. Chemists have worked to find the chemical 
constitution of the vitamin itself and so isolate 
the vitamin from the food in which it was found. 
The physicist has used different methods among 
which the colorimetric and spectrographic are 
the outstanding. 

The purpose of this research was to investigate 
the possible application of the magneto-optic 
apparatus of Allison in the determination of the 
vitamin content of various substances. The ap- 
paratus has been shown to be very sensitive and 
accurate as well as very rapid in making chemical 
analyses. 


EXPERIMENTAL 


The apparatus used for this work consisted of 
a modified form of the magneto-optic apparatus 
used by Allison' in his work and is shown in 
Fig. 1. As a source of illumination, a mercury arc 
was used in connection with a good monochrom- 
ator adjusted to give the 4358A line only. 


1F. Allison and E. J. Murphy, J. Am. Chem. Soc. 52, 
3796 (1930). 


A quartz, potassium-hydride, photoelectric 
cell was used to record the light minima. The 
current from the photoelectric cell was amplified 
by the use of a General Electric FP-54 Pliotron 
tube in a modified DuBridge? circuit. The varia- 
tions in current in the photoelectric cell circuit 
were recorded by a Leeds and Northrup high 
sensitivity galvanometer. The photoelectric cell 
and amplifier circuit including batteries were all 
placed in a closed metal-covered box to shield 
them from stray electric fields. The galvanom- 
eter which was placed directly on top of this box 
was also carefully shielded electrically. 

Only half of the bilateral trolley system used 
by Allison was employed in this experiment. The 
two nicol prisms were held parallel and the mag- 
netic fields in the helices assisted each other; i.e., 
the fields in the two coils both rotated the plane 
of polarization in the same direction. 

The cylindrical, glass cells used to hold the 
liquids in this experiment were 16.3 cm long and 
2.2 cm in diameter. The ends of the cylinders were 
flared and then ground. Glass windows were 
fastened to these flared ends by means of sodium 
silicate. These cells were placed inside the helices 
and were held co-axial with them during the 
experiment. 

After the apparatus was set up, it was cali- 
brated and adjusted for its zero position. This 
was done by using carbon bisulfide in both cells. 
As a check on the apparatus, minima were also 
read with different inorganic compounds in one 


?L. A. DuBridge, Phys. Rev. 37, 392 (1931). 
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Fic. 1. The plate shows the arrangement of the apparatus. The helices, lenses and nicol prisms 
are all mounted on a photometer bench. The movable helix was mounted directly on the photom- 
eter holder. This helix was moved back and forth along the light path scale by means of a long 
screw, which is shown on the plate. The large metal covered box, enclosing all of the photoelectric 
cell and amplifier circuit, is shown at the extreme left. The shielded galvanometer, mounted di- 
rectly on top of this box, is barely visible because the shield surrounding it was painted a dull 
black. Part of the control panel for the photoelectric cell and amplifier circuit is shown on the 
front of the metal covered box. A small reversible, electric motor used to operate the movable 
trolley may be seen in the lower left-hand corner of the plate. The trolleys and electrical path 
scale, not shown in this plate, are immediately above this apparatus and are parallel to the pho- 


tometer bench. 


of the cells. Some of the compounds used in- 
cluded hydrochloric acid, nitric acid, zinc 
chloride, barium chloride, water, ferrous sulphate 
and ferric chloride. These solutions were all very 
dilute. The minima found for all these compounds 
checked very closely with those found by Allison. 

In order to find whether there were distinct 
minima for vitamin A, one of the cells was filled 
with pure ethyl alcohol and a drop of Haliver oil 
was added. (The Haliver oil, which had 100 times 
the vitamin A potency of cod-liver oil, was ob- 
tained from Park, Davis and Company of De- 
troit, Michigan.) Minima were read on this solu- 
tion as the movable trolley was moved from 0 to 
41 along the electrical path scale; i.e., each time 
the galvanometer showed a deflection, the scale 
reading was recorded. This procedure was fol- 


lowed five times and the scale readings were then 
rechecked. Only those scale readings were con- 
sidered which could be rechecked consistently. 
Then the tube which contained the Haliver oil in 
alcohol was thoroughly cleaned and filled with 
pure alcohol alone. Five sets of minima readings 
were then taken and rechecked as before. The 
two final sets of scale readings were then com- 
pared and it was found that there were several 
scale readings which were characteristic of the 
Haliver oil alone. 

Cod-liver oil was next tested in the same man- 
ner as Haliver oil and minima were read which 
were characteristic of this oil. In like manner, 
other substances containing vitamin A were tested 
and it was found that all of the substances which 
were known to be rich in vitamin A content gave 
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MAGNETO-OPTIC 


a minimum reading at 32.9 on the electrical path 
scale. This minimum was therefore concluded to 
be related to vitamin A. 

Since irradiation of cod-liver oil with ultra- 
violet light reduces the vitamin A potency of the 
oil, a sample of the cod-liver oil used was ir- 
radiated under a quartz mercury vapor lamp for 
four and one-half hours at a distance of one foot. 
This sample was tested in the apparatus and it 
was found that the minimum at 32.9 showed up 
only very faintly. 

A sample of oxidized cod-liver oil was obtained 
by bubbling air through it for fifteen hours at a 
temperature of 92+3°C. When this sample was 
tested, it showed no minimum at 32.9. 

Peanut oil and Wesson oil were also tested, 
but neither of these two showed a minimum at 
32.9. Different carotene solutions were also 
tested. Carotene crystals and a carotene solution, 
both of which were obtained from Mead Johnson 
and Co., did not show minima at 32.9. One sample 
of carotene dissolved in Wesson oil, which had 
been obtained from the U. S. Department of 
Agriculture, did show a minimum at 32.9; but 
this was probably due to slight traces of vita- 
min A. 

Some carotene crystals were dissolved in alco- 
hol and the solution was irradiated under a 
quartz mercury vapor lamp for a period of four 
and one-half hours at a distance of one foot. 
When this irradiated solution was tested, it 
showed the minimum at 32.9. A similar solution 
which had not been irradiated did not show this 
minimum. 

RESULTS 


A list of the substances tested and the solvent 
used is shown in Table I. This table also shows 
whether the minimum at 32.9 appeared for these 
various substances. 


DISCUSSION 


Since the minimum at 32.9 was absent in sub- 
stances which did not contain vitamin A and 
present in those substances which did, it was 
concluded that this minimum is closely related 
to, or associated with, vitamin A. The fact that 
the minimum was weakened or destroyed by 
irradiating or oxidizing the cod-liver oil gave 
additional strong evidence that it was closely 
related to the vitamin A content of the oil. Other 
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TABLE I. 
Minimum 
Substances Solvent at 32.9 

— alcohol no 
water no 
Haliver oil alcohol yes 
cod-liver oil alcohol yes 
oxidized cod-liver oil alcohol no 


irradiated cod-liver oil alcohol very faint 
irradiated cod-liver oil water very faint 
spinach juice alcohol yes 
spinach juice water yes 
orange juice alcohol yes 
orange juice water yes 
tomato juice alcohol yes 
tomato juice water yes 
peanut oil alcohol no 
peanut oil water no 
Wesson oil alcohol no F 
Wesson oil water no 

egg yolk water yes 
carotene (1) alcohol no 
carotene (1) water no 
carotene crystals (2) alcohol no 
irradiated carotene crystals (2) alcohol yes 
carotene in Wesson oil (3) alcohol yes 


Carotene (1) solution was obtained from Mead, Johnson 
and Company, Evansville, Indiana. It had 10 times the 
antixerophthalmic potency of ordinary cod-liver oil. 

Carotene (2) crystals were obtained from Mead, Johnson 
and Company, Evansville, Indiana. 

The carotene (3) in Wesson oil was obtained from the 
U.S. Department of Agriculture, Bureau of Chemistry and 
Soils. 

The irradiated carotene crystals (2) were irradiated under 
a quartz mercury lamp for four and one-half hours at a 
distance of twelve inches. 

The egg yolk solution was obtained from the Poultry 
Husbandry Department of lowa State College and it was 
known to contain vitamin A. 

The irradiated cod-liver oil was obtained by irradiating 
Norwegian cod-liver oi! under a quartz mercury vapor lamp 
for four and one-half hours at a distance of twelve inches. 

The spinach juice used was obtained from a can of Del 
Monte spinach. 

The orange juice and tomato juice were obtained directly 
from the fresh ripe fruit. 


minima were also examined, but the minimum at 
32.9 was the only one for which the evidence was 
sufficient to indicate a close correlation with 
vitamin A. 

This minimum was not found in any of the 
carotene solutions except the one carotene- 
Wesson oil solution. The pure carotene crystals 
(2) and carotene solution (1) showed no minimum 
at that point. 

Much work has been done in recent years to 
find the relationship of carotene to vitamin A. 
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Moore*: ‘found that carotene behaves in vivo as a 
precursor of vitamin A. von Euler, von Euler, and 
Hellstrom® used the antimony trichloride test on 
carotenoids and found that they gave a reaction 
similar to that due to vitamin A, and that caro- 
tene from carrots acted the same upon the growth 
of their test animals as vitamin A. 

Bowden and Snow° irradiated a solution of 
carotene in cyclohexane in an atmosphere of 
nitrogen. Their exciting light was the 2650A 
line of mercury. After a few hours of irradiation, 
they found that the solution had a strong ab- 
sorption band at 3280A, which is characteristic 
of vitamin A materials. The solution also gave 
the blue coloration when treated with antimony 
trichloride in chloroform. From these tests, they 
concluded that carotene was changed to vitamin 
A by the action of ultraviolet light. However, 
their conclusions were questioned by Heilbron 
and Morton.’ 

In the experiments described here, carotene in 
an alcohol solution was irradiated with ultra- 
violet light; and it was found that, after irradia- 
tion, this solution gave the minimum at 32.9; 
whereas, before irradiation, this minimum was 
not present, indicating that a small amount, at 
least, of carotene can be transformed into vitamin 
A by irradiation. 


*T. Moore, Biochem. J. 23, 803 (1929), 

*T. Moore, Biochem. J. 24, 692 (1930). 

5 Beth von Euler, H. von Euler and H. Hellstrom, Bio- 
chem. Zeits. 203, 370 (1928). 

®F, P. Bowden and C. P. Snow, Nature 129, 720 (1932). 

*I. M. Heilbron and R. A. Morton, Nature 129, 866 
(1932). 
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A sample of egg yolk solution was obtained 
from the Poultry Husbandry Department and 
tested in the apparatus. This sample was part of 
the solution used in feeding experiments and it 
was known to contain vitamin A. This sample 
also showed the minimum at 32.9 when tested in 
the apparatus, thus checking the spectroscopic 
investigations by von Euler and Klussman* who 
found the absorption bands in egg yolk solutions 
characteristic of vitamin A. 


CONCLUSIONS 


Conclusions which may be drawn from the 
results of this investigation are as follows: 


(1) The magneto-optic apparatus can be used to de- 
termine the vitamin content of various substances, 

(2) The photoelectric cell may be used successfully 
with the magneto-optic apparatus. 

(3) The minimum at 32.9 on the electrical path scale 
used shows evidence of being very closely related 
to the vitamin A content of different substances. 

(4) Pure carotene in an alcohol or water solution does 
not give the minimum at 32.9, but after irradiation 
with ultraviolet light, it does give this minimum. 

(5) The magneto-optic method is very rapid in deter- 
mining the vitamin content of different substances. 
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Mutual Impedance of Grounded Wires Lying on the Surface of the Earth when the 
Conductivity Varies Exponentially with Depth 


Marion C. Gray, American Telephone and Telegraph Company 
(Received May 18, 1933) 


This paper presents a formula for the mutual impedance of any insulated wires of negligible 
diameter lying on the surface of the earth and grounded at their end-points, on the assumption 
that the conductivity of the earth varies exponentially with depth. Various special cases are 
briefly discussed. 


FORMULA established by Foster! for the mutual impedance of grounded thin wires lying on 
the surface of a uniformly conducting earth has now been extended to the case of an earth 
whose conductivity varies exponentially with depth. The formula obtained may also be regarded as a 
generalization to finite wires on the surface of the earth of a mutual impedance formula recently 
published? for infinite wires at or above the surface of an earth of exponentially varying conductivity. 
As in that paper, we assume that the conductivity of the earth at depth z is defined by the formula 


y=ye*, 230, (1) 


where yo is the conductivity at the surface (the x—y plane). 

The further assumptions made are, as usual, that the conductivity of the air is zero, that the 
capacitivities of the earth and of the air are both zero and their inductivities both equal to that of 
free space v, where v has the value 1.256 X10~* henries per meter. Finally, all displacement currents 
have been neglected, so that the frequency w/2z is restricted to relatively low values. 

With these assumptions the differential equations satisfied by the electric field in the earth are 


=T (2) 
(3) 
(4) 


where I) = (iwvyo)! is the propagation constant at the surface. For positive values of b a solution of 
Eq. (4) which vanishes at an infinite depth is* 


ff Pw vje*?K sin ux cos vy du dy, (5) 


where 


and F, is an arbitrary function of its arguments. 
Substituting this value for EZ, in Eqs. (2) and (3) the following solutions for E, and E, are obtained: 


1R. M. Foster, Bull. Am. Math. Soc. 36, 367-368 (1930); Bell System Tech. J. 10, 408-419 (1931). 

*M. C. Gray, Physics 4, 76-80 (1933). 

§ For the Bessel functions J,(z), K,(z), J,(z) we follow the definitions given by G. N. Watson, Theory of Bessel Functions, 
pp. 40, 78, 77, Cambridge, 1922. 


mw 2uBe>=!2 
Bam Ka(s) | c08 ux cos vy dude, (6) 
ovo 
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2vBe>=!2 
ff Fyi(u, Ap, Kasilf)+— sin wx sin vy du dp, (7) 
1 


where 7 = 28(u?+»’)!, and the functions F, and F, are again arbitrary. For negative values of b 
similar solutions are obtained, with the K functions replaced by the corresponding J functions. 

Starting from the solutions (5), (6) and (7), in place of the solutions of the wave equation used by 
Foster, the impedance formula is obtained by the second method outlined in his paper, and it is found 
that the mutual impedance of wires S and s lying on the surface of the earth and grounded at their 
end-points may be expressed in the usual form 


Zu= f /dS ds+cos eN(r) WS ds, (8) 


where the integration extends over the two wires S and s;r is the distance between the elements dS 
and ds, and ¢ is the angle between their directions. The functions Q(r) and N(r) are defined as follows: 


O(r) = (1/2070) F(u) Jolru)du, (9) 
= (iwv/2aT) de. (10) 


where 


F(u) = 0), O58, 
= I_29,(— — 261 0), b=0, 
F(u) = (To/#) CK a41(26T 0)/Ka(281 —f(u)J—2Bu/(a—1), O05), 
a= 


In these formulae the distances r, S, s are measured in meters, the impedance Z,2 in ohms and the 
conductivity 7 in mhos per meter. 4 
The following special cases of the emai may be noted. 
(1) If one wire S is straight and infinitely long, the mutual impedance gradient parallel to S and 
at a positive distance y from S is 


dZ32/dS = Gin #00) COS Vu du, 0=), 
0 


= 281 o) 261») ] COS Vu du, b=0, 
0 


a formula which has already been established.* 
(2) If bis allowed to become negatively infinite while at the same time 7») becomes infinite in such 
a way that y:= —70/) remains finite, Eqs. (9) and (10) become 


Q(r) = (1/471) [Ho(2iwryir) — Yo(Ziwvyir) J, 
N(r) — ((twv)? 71/16) diwvyir) — Yo(Ziwvyir) |, 


* Reference 2, page 77. 
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where Hy(z) is Struve’s function of zero order, and Yo(z) is the Bessel function of the second kind of 

zero order as defined by Watson.® These values agree with the mutual impedance formula for an earth 

consisting of a conducting layer at the surface only, as already obtained by Riordan and Sunde.® 
(3) For small values of 6 the first terms in the expansions of Q(r) and N(r) in powers of b are 


Q(r) =1/2ayor — or)Ko(3T or) +: 
N(r) — (1+T + LK o( Por) 07) K of) — 07) K o( 31 or) 


The terms independent of } in these expansions give precisely Foster’s formula for a uniformly con- 
ducting earth.’ 


(4) If I) is small compared with 6 the first terms in the expansions of Q(r) and N(r) in powers of 
w are 
Q(r) = (1/2270) [(1/r)e7'? — 30Ei (37d) J+ (iwvbr?/322) log w+0(w), 0<b, 


= — —c log r]+0(w), b=-—c<0, 
N(r) = (iwr/4rr)(1 — (iwvryo/2b) log w+0(w) 


where Ei(z) is the exponential integral as defined by Whittaker and Watson. 

In order to obtain convergent integrals in the expansion for Q(r), terms independent of r have 
been added to the integrand and, after integration, all such terms have been neglected, since they 
contribute nothing to the mutual impedance. The first term in the expansion gives the direct-current 
mutual resistance, which has been verified independently.® 

The term in w log w, occurring in the expansion of Q(r) for positive values of b, is worthy of men- 
tion as it appears to be new. In all previously published cases the expansion in powers of w has been 
found in the form 


where R is the direct-current mutual resistance and N the direct-current mutual inductance of the 
wires. Reference should be made, however, to the anomalous case mentioned in Riordan and Sunde’s 
paper’? on the two-layer earth, when the lower layer has an infinite conductivity, a condition which 
has a formal resemblance to the exponential case when 6 is large and positive. It is then found that 
two distinct values of the mutual inductance WN are obtained, depending on the order of proceeding to 
the limits w=0, A»= «. A fuller discussion of this point will be found in the paper cited. 


5 Reference 3, pages 328, 64. 

6 J. Riordan and E. D. Sunde, Phys. Rev. 37, 1369-1370 (1931); Bell System Tech. J. 12, 162-177 (1933). 

7 Reference 1, page 411. ; 

SE. T. Whittaker and G. N. Watson, Modern Analysis, Third Edition, p. 352, Cambridge, 1920. 

® This formula for the direct-current mutual resistance for positive values of 6 was obtained by us in 1930 for use in 
connection with test measurements of the earth's conductivity, but it was not published at that time. 

10 Reference 6, pages 165-166. 
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